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ilBSTRACT 


lhl«  papor  presents  s  new  nunerlosl  method,  the  error  method, 
for  solving  psxsb<aio  type  pertlsl  differential  equations,  linear  or  non¬ 
linear.  In  particular,  ty  conparlng  with  the  regular  successive  iterative 
method,  more  lieneflolal  results  In  the  applloatlon  to  non-linear  problems. 
Three  non-linear  examples  were  studied  by  using  this  method.  All  resulted 
In  large  reduction  of  number  of  Iteration  loops  and  CFV  time  required  In 
comparing  to  the  corresponding  regular  suocesalve  method  used. 

Generalisation  and  modification  of  this  method  appears  appropriate 
to  extend  Its  application  to  elliptical  type  partial  differential  equation 
so  that  problems  with  "Isolated"  events  (  such  as  those  with  Ignition  spots) 
nay  be  handled  with  this  method. 
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I.  Introduction 

Ihla  paper  preaenta  a  new  mMerioal  aethod,  the  error  aethod, 
ror  aolYlng  paiabollo  type  partial  differential  e^uatlona,  linear  or  non¬ 
linear.  Obaervinc  froa  the  oaaea  atudled,  it  appeara  that  the  application 
of  the  error  aethod  to  non-linear  proUeaa  reaulta  in  aore  beneficial  than 
to  linear  probleaa  in  ooaparing  with  the  conventional  workable  aucoeaalve 
iterative  nunerloal  achene.  Ihia  new  method  featurea  wide  applicability. 
Ita  coaputational  achene  oonvergea  faat.  The  following  aeetlona  will 
deacrlbe  the  fundanental  derivation  of  thla  method,  ita  applioations  to 
non-linear  probleaa  and  the  future  potential  developnenta. 

II.  The  Suooeaaive  Iterative  Hethod 

To  conatruot  the  baalc  numerioal  achene  of  the  error  method, 
let'a  etart  by  reviewing  the  conventional  auooeaaive  Iterative  method. 
Conalder  the  problem  of  am  one-dineaaional,  iinateady  heat  conduction  in  a 
aolid  plane  wall  with  finite  thiokaeaa.  Ita  governing  eqmtion  ia, 

Tt  -  Txx  .  (1) 

where  T,  the  uakaoMa  temperature  field  and  x,t,  the  effective  apacioua 
and  tine  coordiaatea.  The  aet  of  initial  and  boundary  condltlona  may  be, 

t  -  0,  T  -  T^(x) 

X  •  0,  T  •  C, 

X  -  1,  T  -  Ca  , 


Mher*  1,  the  effective  well  thiokBees  aad  C| ,  Ci,  two  fixed  teepexetuxe 
values  os  the  two  wall  suxfaces.  Tq(x)  xepxeaeata  the  initial  tespexsture 
profile  aoross  the  wall  thiokneaa. 


To  such  a  prohles,  the  finite  difference  suooessive  iterative 
nunerloal  method  oan  be  used.  Ihe  field,  along  x  direction,  is  first 
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Figure  1.  Grid  Systes  for  the  Nxinerioal  Scheme 


divided  into  (N-1)  grids.  Each  grid  sises  ax.  An  implicit  finite  difference 
equation  is  constructed  from  equation  (1).  i.e. 


—  •  ■ 


a-t 


+ 


(3) 


where  FtT  )  _ 


and  n  varying  from  2  to  N-1.  The  solution  of  equation  (3).  together  with  the 
two  boualary  ooaditloas  specified,  determine  the  vmrioua  T^,  ^as  are 


(3) 


giv«n«  Ulth  sueocssiv*  itezatlv*  nethod,  solutloa  of  •qmtlon  (3)  starts 
tgr  aasuaing  a  aet  of  values »  together  with  the  specified  iaforaation 
for  F(Tn^»)  froa  tioe  iastanttS)  auhetitutiag  into  the  right  side  of 
equation  (3)»  then  calculating  a  new  set  of  Trvu^i,,^ .  'Ihis  new  set  of 
is  then  xeplaolng  the  old  set  of  to  plug  into  the  ri^t  side  of 

equation  (3)  again  and  calculating  the  next  new  set  of  .  Such 

procedure  repeats  itself  until  the  calculated  set  of  revealing  no 

signifioant  difference  froa  the  previous  set  of  .  Ibe  last  set  of 

calculated  values  are  considered  the  satisfied  solution.  Ihe 

selection  of  the  starting  profile  is  not  unique.  Ordinarily,  to  use 

the  profile.  .  as  the  first  trial  profile.  appear  not  a  had 

choice.  Since,  intuitively,  as  we  nay  expect  that  the  correct  unknown  profile. 

I  should  not  be  different  froa  the  known  profile  of  the  previous  tiae 
instant  .Tm.ht  »  too  much  if  ^t  is  saall  enough  and  no  ’’isolated”  event 
occurs.  To  assure  the  stable  convergence  of  the  calculation  steps  and  lead 
to  reasonable  solution,  it  appears  that  the  criterlan. 

A-t 


ought  to  be  observed  in  selecting  the  relative  sises  of  and  ^t. 

According  to  above  description,  a  Fortran  program  (SPDl).  atta^ed 
as  Appendix  A.  was  eonatruoted  to  allow  digital  computer  to  perfora  the  steps 
la  the  iterative  coaputatioa  prooed\ire  for  this  one  dlneaslonal.  unsteady 
heat  conduction  problen,  A  oase  with  the  initial  condition. Tg(x}.  poasesaing 


oa»  surfao*  ivaperatui*  alevatad  level  ms  studied,  Soae  oslculated 
tempera ture  profiles  at  various  tise  instants  axe  shewn  in  IhliLe  1.  In 
this  table,  the  auaber  of  iteration  loop  required  for  calculating  each 
tenpexature  profile  is  also  indioatedi  so  is  the  total  CPU  tine  used. 
Also,  oases  with  highly  non-linear  initial  conditions  were  also  investi¬ 
gated  with  this  program,  the  results  indicate  that  the  successive  method 
appears  very  effective  to  handle  such  linear  problems. 


Table  1  Some  Tenpexature  Profiles  for  a  Qne-Dimensioaal, 

Unsteady  Heat  Conduction  Problem  (Successive  Method) 

Initial  Temperature  Profilei  7^00,0  at  x-0,0  and  t-0,0 

k0,0  at  X  0.0  and  t"0,0 
Boundary  Conditions!  T^OO.O  at  x"O.0  and  t  0.0 

T-  40.0  at  x-1.0  and  t  0.0 

The  label!  Orderly  number  of  Temp  Profile,  Ho.  of  Iteration  Loop,  Time  Inst, 
Values  of  Temperature  Profile 


1 

0  0. 

00000 

200.00000 

40.00000 

40 . 00000 

40.00000 

40.00000 

40.00000 

40.00000 

40 . 00000 

40.00000 

40.00000 

40.00000 

40.00000 

40.00000 

40 . 00000 

40.00000 

40.00000 

40.00000 

40.00000 

40.00000 

40.00000 

40.00000 

20 

8  0. 

01900 

200.00000 

167.54154 

137.17704 

110.59675 

88.81256 

72.08412 

60.03V»0l^ 

51.88440 

46.69865 

4  3. '590 19 

41.83147 

40.89040 

40.41.3.1.4 

40.18326 

40.07784 

40.03172 

40.01242 

40.00468 

40.00168 

40 . 00053 

40.00000 

40 

9  0. 

03900 

200.00000 

177.23933 

155.20179 

134.54165 

115.78639 

99.29813 

Bfj .  ;;>rj89o 

73.67917 

64 . 42'334 

57.25884 

51.87905 

47.96330 

4f..l9H94 

43.30546 

42.04644 

41.23314 

40.72176 

40.40711 

40.?14'36 

40.09160 

40.00000 

60 

9  0. 

05900 

200.00000 

181.46837 

163,32715 

145.94209 

129.63219 

114.65246 

101 . 18298 

89.32492 

79.10333 

70.47551 

63,34362 

57.56936 

'32.98915 

49.42818 

46.71205 

44.67555 

43.16817 

42.05653 

41.22417 

40.56913 

40.00000 
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III  The  Error  Method 

In  applying  conveaslMal  successive  nuserical  scheme  for  solving 
parabolic  partial  differential  e<iuations,  ve  may  •‘feel**  one  feature,  l.e,, 
the  "passlveness”  of  the  oonputatlon  scheme.  Except  selecting  the  slses  of 
the  Independent  variables  and  specifying  the  initial  condition,  there  is  no 
more  handle  which  can  be  used  to  alter  or  control  the  development  of  the 
computation.  Fuxrthermore,  the  convergent  rates  of  the  successive  method 
axe  often  small,  especially,  to  non-linear  problems.  Frequently,  hundreds, 
or  even,  thousands  iteration  loops  are  required  to  achieve  a  single  step 
solution.  On  account  of  all  these  unsatisfactory  characteristics  with  the 
coQvensional  successive  method,  the  error  method  is  developed.  Essentially, 
the  major  difference  between  the  successive  and  the  error  method  is  that  the 
later  imposes  mechanism  for  estimating  and  correcting  "errors"  of  trial 
values  of  the  unknown  quantities  in  a  regular  successive  routing.  So  the 
overall  computation  procedure  is  greatly  shortened.  To  illustrate  the  basic 
scheme  of  the  error  method,  let  us  re-conslder  the  unsteady,  one-dimensional 
heat  conduction  problem  of  last  section.  Let  ^  represents  the  first 
trial  temperature  profile  at  m+l  tine  instant,  represents  the  "true 

solution"  (  still  unknown  yet).  It  is  clear, 

Ttt.,.  -  (5) 

where  <5^,  the  devlatim  of  trial  profile,  ,  from  the  "true  solution", 

•  To  start  the  iteratitti  process  by  substituting  T^'',n  into  the  right 
side  of  Equation  (3) t Just  as  same  as  done  in  successive  method,  and  decomposing 


(6) 


aooordiag  to  equaUon  (5),  we  obtala 


<1 


,  —  '  -T-  >  S'v,-j-l'+  ^n-l 

+  ’'^^-,*^-1^7^  .  {«) 


where  ^  ,  the  oalculated  tenpexatuze  profile  at  «+l  time  Instant  based 
on  „  .  Substractlng  both  sides  of  equatiwi  (6)  by  „  ,  we  then  obtain. 


TCl  > 


^  _  ■  Srn-i-t- 


.  (7) 


Equation  (?)»  the  error  equation,  provides  N-2  required  equations  for 
solving  set  of  N  unknowns,  (Ts.  The  other  two  equations  required  are  fomed 
from  the  two  boundary  conditions.  For  this  particular  case,  it  is  obvious 
that  determining  the  set  of  S'  values  fro*  equation  (?),  the 

^loes,  ,  oan  then  be  Modified  and  iteration  processes  repeat, 

if  needed,  until  the  aatlafactory  "true  solution"  is  achieved.  To  solve 
equation  (7)  »na  tne  cor»sponding  oounoary  conditions  for  S',  a  set  of 
reeurreney  forauiar  are  derived.  Let 


S|ni"  Bm  , 

where  and  are  undetermined  correlation  parameters  between  neighboring 
Re-arrange  equation  {?)»  we  obtain. 


2  (  14 


(7) 


KllBlMtlng  from  a^tmtlon  (9;  and  equation  (8;,  we  hawe  the  leeorreney 
ronulart 


•n+i 


,  »  I 


B 


(10) 


To  »xmn  the  opexatloa  oa  the  zeouzxeney  foxaular,  we  kaew.froa  equation  (9; 
and  boundary  eondltlon,  U,  that 


uu 

Bs  -  -2(  A+  ^ 

After  all  the  A^,  and  B^  ,  for  3$n;sN,  are  deteralned,  since  S"  0, 
therefore « 


Ihe  rest  ^a  axe  oaloulated  froa  equation  (8)  In  reverse  order.  The  deteralned 
<r values  axe  used  to  aodlfy  the  trial  values  of  unknoan  variables,  , 

and  the  Iteration  cycles  repeat,  if  needed,  until  satisfactory  "true  solution” 
Is  obtained.  It  is  clear,  to  this  problem,  a  linear  case,  equation  (?)  la  a 
linear  equatlM  of  .  However,  for  non-linear  problea,  equation  (7)  would 
be  a  non-linear  equation  of  ^  .  To  facilitate  the  calculation,  linearisation 
procedure  Is  used  to  convert  equation  (7)  Into  linear  fora  of  The  other 
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8t«pB  Of  derivation  are  Identical  to  those  described  above. 

For  clarity,  the  above  description  about  the  derivation  of  the 
basic  error  aethod  will  be  suanarlsed  as  follows  t 

(1)  Corresponding  to  the  governing  differential  e^untloas  as  well  as 
the  boundary  conditions,  fom  the  inplloit  finite  differenoe  equations, 
ready  to  start  the  regular  iteration  procedure. 

(2)  Sy  using  the  first  trial  values  of  the  unknown  variables,  say 
caloulate  the  unknown  variables,  say  . 

(3)  By  reoognielng  .idiere  .  the  "true  solution", 

deduce  the  linear  error  equatlons.Ry,*  f(  )f  where  R„-  # 

fron  the  finite  difference  niwerlcal  scheme. 

(4)  Let  .  where  .  B„*,  ,  the  undetermined 

correlation  parameters.  Derive  the  recurrency  formular  for  .  B„^^  , 

l.e. . 

A  M+i  “  G(Aki  ) 

Bw+t  “  H(A«  .  B|^  ), 

from  the  error  equations.  Determine  the  starting  A  .  B  values  from  the  related 
boundary  condition. 

(5)  Calculate  A„  .  B.,  for  n  up  to  N.  Determine  first,  with  the  help 
of  <Mie  boundary  condition  and  then  with  the  correlation  fomular. 

(6)  Modify  and  repeat  steps  (2)  to  (5),  if  needed,  until  the 

satisfactory  "true  solution"  is  obtained, 

VI th  above  procedure,  a  Fortran  program. SPD2.  was  constructed  and 
attached  as  Appendix  B.  Oie  case,  sane  as  studied  by  SFDl.  was  investigated 
by  this  program.  The  results  are  identical  to  those  fron  SFDl.  Thble  2  shows 


(9) 


the  eoaparisoQ  of  the  number  of  itermtlon  loops  CPU  time  required  by 
SP2)2  and  SPDl  for  that  particular  case  studied.  From  Tbble  2,  It  appears 


Ifcble  2.  A  coaparlson  of  the  Number  of  Iteration  Loop  and  CPU  Time 
Required  by  Successive  and  Error  Methods  for  the  Problem 


Studied  In  Tkble  1 

No. of  iteration  loop 

Total  CPU  tine  required  for 

Program 

per  time  step 

So  time  steps 

SPDl 

6—9  (aoot,9) 

5.30  ssc 

SPn2 

1  (occasionally  2, 

^.53  sac 

once  3) 

that  the  CPU  time  saved  by  the  error  method  In  reference  to  successive 
method  Is  only  about  20  %,  ».  value  mot  very  Impressive.  HoMever*  this 
observation  may  also  be  Interpreted  as  that  the  successive  method  Is  already 
fairly  effective  Im  applying  to  linear  problem. 

IV  Ihe  Application  of  Error  Method  to  Non-linear  Problems 

Three  non-linear  exanples  are  presented  to  Illustrate  the  appllca- 
tlcm  and  effectiveness  of  the  error  method.  To  facilitate  the  presentation, 
brief  qualitative  description  of  these  three  problens  are  given  first  as 

fOllOHSI 

(1)  A  steady  laminar  tubular  flow  heat  transfer  problem  with  fluid 
themal  conductivity  varying  linearly  with  local  temperature.  This  problem 
features  with  non-linear  heat  conduction  term  In  the  transveroe  dlreotlcm 
to  the  tube  flow. 


Figur*  2.  Ihe  Coaxlftl  Tubular  Flow  Sketch 


Initial  velocity  profile  is  discontinuous  at  inner  tube  wall  location.  The 
omeentxatlon  profiles  of  both  the  inner  and  outer  tube  gases  are  disconti¬ 
nuous  at  the  inner  tube  wall  location.  All  the  difusslon  terns  in  the 
transverse  direotion  are  non-linear. 

(3)  An  unsteady «  oae-dinenslonal  fin  problen  with  radiation.  The 


Figure  3*  Ihe  Sketch  of  a 
radiative  fin 


scheaatic  dlagxan  is  shown  in  Fig\irs  3.  To  this  problen,  the  highly  n<m- 
linear  radiation  tern  presents  in  both  differential  equation  as  well  as  the 
boundary  condition. 


Th«  foraulatlm  of  all  three  probleas  described  are  shorn)  in 

Ihbls  3. 


Thble  3.  Forsulatlon  of  the  three  Non-linear  Frobleas 


Items 

Tube  H.T*  Problem 

Coaxlle  Tube  Problem 

Hadiatlve  fin  I^oblem 

Governing 

Equations 

1 

rr-i 

e.X:^ 

Dh  (-r  C,C©^4Tr©V6Tr  ©t^V) 

c,3  ed'CT.'T.o^y-^ 

%.= 

Initial 

and 

Boundary 

Cttidltions 

1 

1 

X^J  -T-To(.ry 
=  o 

It -0.5,  T  ~Tw 

T=TaC'r)Y»V'r^/T4-V’^ 

'  *  i'T  7?^  »r 

-t-o^  9--! 
e=-i 

C,  =  -^~  fcioV  v**Jiie^ 

<•  -mwx  ‘  _,, 

Finite 

Difference 

equations 

j 

,,  S-^X^(isL?r 

‘^a*'=Ct-1nA--0 

Tv^-Tvv 

Temperature  expressions 
are  identical  to  those 
in  left  column • 

+^<V+V> 

J.  y 

C,J  '•<) 

ffx»r*-RaSc 

C3‘'i>/CT-l> 

'V  c  <v-  ^  Yo  r>  Yf. 

e  .  ®"'  «  ^  e  . 

<^4,  *  5^x^-»-  i-  C3  J>,i  - 

C  * -f- ~ -f- r  ■^\ 

"HJe  ^<4X  •(tk—  N 

(Continue  at  next  F*^ge) 


(12) 


3.  (contiauad) 


liana 

Tuba  H.T.Froblan 

Coaxila  Tuba  Problan 

Rxdiativa  Fin  Problan 

Linanriaad 

Error 

Equations 

ranpazatuxa  arror  aqa, 
bra  Idantical  to  thoaa 

La  laft  colunn, 

c- 

^  C  ' 

v/v 

Su’  O.XS  Cc-gr  • 

^  ^30i^c«w-ii  ^  ft.'  ^  0.i.MH:Cei^ 

^  ‘1>  lb 

Cp5  =  0,(1  ®.\;  ^  V©,-  +jV’ 

Bacurranoy 

foxnular 

e.  =i?i.cR. 

V'  A; 

^,=  ^ 

<AX 

T  <v^X  Y»  (^eSSis. 

vl;-, 

Infomntioa 

darlvad 

fron 

boundary 

oondltlons 

s  =-|i 

H-i  Am 

O  Sm 

t-rr  Yi*  t 

Ax=l-».  fri-o.a. 

Y*  ^  . 

'K  0  _  £i££. 

J  -lli - hL 

-  A 

(eontlnuad  on  naxi  pnga) 


(13) 


TtelAe  3«  (oontlBusd) 


Tube  H.T.Froblen 

Coaxlle  Tube  Problen 

Radiative  Pin  Problen 

Modified 
Ttlal  Profll 
for  next 
IteratlM, 

If  needed 

V  T-j-r 

Y') 

Nonenolatun 

unspec* 

V 

above  aztd 
left  expre* 
sslons, 

?,fluld  density 
u, fluid  velocity 

I, fluid  tenpexn- 
ture 

Xytube  axlle  eoor, 
r,tube  rad  ooor. 
Re, Reynolds  nuabez 
Pr.Fxandtl  nunber 

tenperatuce 
Profile 
tuke  wall 
teamzatxire 
j  (or  J),jth  node 
point  along  r* 
direction 

N, total  nunber  of 
nodal  points 
alMg  r-diree* 
tlon 

deviation  of 
trial  value  froi 
"true  solution.’ 

Yg, outer  tuke  gas  cone 
1^, Inner  tube  gas  cono 
Y^(r),  entrance  Y,, 
profile 

Yf.^(r),  entrance  Y^ 
Profile 

Sc,  Schnldt  number 

.Tj^, Initial  and  fin 
•  root  tenperature 

Tm  ,environnent  tenp. 

£  ,ennlalvlty 
d*  ,Stef-BoltBnann  Cons, 
h,  eOBv.  h.  t.  Coef. 

1,  fin  length 

P,  fin  cross  section 
perimeter  length 
k,  fin  themal  oond. 

A,  fin  cross  section 
area 

(14) 


BRsed  upon  fontular  In  Ikble  3*  oonputer  progsann  with  error 
Method  were  designed  for  the  above  three  probless  as  SPD6,SPD4  and  SPD6. 
Ihey  axe  attached  as  Appendix  F,D  and  H.  Ihree  corresponding  prograas  ( 
SF05»SFD3  and  SPD7  )tb&sed  on  successive  aethod,  for  these  three  problems 
axe  also  constructed  and  attached  as  Appendix  I,C  and  C.  Conputatlonal 
effectiveness  of  these  programs  were  investigated  on  several  sample  oases. 
lAble  4  shows  several  ooaputational  effectiveness  of  variotis  programs  run. 


Tbble  4  Conparlscm  of  CPU  time  and  Number  of  Iteration  Loop  xeq,uired 
l]gr  Successive  and  Error  Method  for  the  three  non-linear 
problems 


Items 

Tube  Flow 

Coaxlle  Flow 

Radiative  Fin 

Given 

Conditions 

Ro-500.00 

Pr-1.0 

rw-3.0 

H  (total  steps 
oaloulated)" 

30 

Re-500.00 

Pr-1.0 

Sc-1,0 

V3»0 

N-41 

M-40 

nl-1.0 

C  -2.0 

C  -0.2 

T  -0,4 

sffln  th,)-0.05 
N-41 

M-eo 

Program 

used 

Succ 

SPD5 

SPD3 

SPL? 

Brroi 

SPD6 

SPD4 

SPD6 

No,  of  Loops  used 
per  step 

Succ 

49-73 

85-127 

7-9 

Erroz 

3 

4 

4 

Total  CPU 

Tine 

Succ 

53.85  sec 

67,70  sec 

29.05  sec 

Erroz 

5.55  ssc* 

30. 18  sec 

4.71  sec* 

Sprint  only  1/lOth  as  ahny  lines  as  did  by  t 

he  comparable  prograas. 

Tbble  4  indicates  that,  by  using  error  method,  the  CPU  tine  used  can  be 


reduced  to  about  half  or  less  in  reference  to  those  correspoDding  oaleulation 
done  by  successive  aethcd. 


7  DisouMlon  mad  Conoliwlon 

(1)  InTMilgailOB  Is  previous  MctioM  hM  Inliostad  that  th«  arror 
Mthod  devalopad  up  to  this  sts«s  doss  provids  us  ss  sltsxastlva  wsy  fear 
sOlviBs  psssbollo  partial  dlffaraatla.1  aquations,  lha  utilisation  of  this 
nathod.  la  particular  to  aon-linsar  prohlans,  rasuits  la  slgalfloaat  raduelaiT 
CPU  tlaa  and  aunhar  at  Itaratlon  loops  roquirad.  This  faatura  ,  togathar 
with  Its  wida  appilosBliltj  nay  turn  this  asthod  Into  aa  attraetlwa  tool. 

(2^  Tha  raourraaey  taehalqoa  auggastad  hara  oaa  ba  eoasldarad  as  an 
afiaetlva  nathod  tar  solTlag  tha  thraa  Has  aatrix  . 

(3)  Tha  applicability  of  tha  arror  nathod  to  ordlasay  dlffaraatial 
aquations,  llnaar  or  aon-linaar.  la  obvious. 

(4)  To  liaaar  parabolic  typa  aqisitlon,  the  above  reourxeaey  taonaique 
can  ba  applied  directly  to  the  tmknowQ  parameters  lastaad  or  the  errors  oi 
the  trial  values  oif  the  unknown  pazaaetara.  It  thus  provides  another  alter¬ 
native  method.  To  non-linear  equation,  however,  the  direct  application  of 
tha  racurrency  technique  to  the  unknown  parameters  auiy  not  ba  abla  to  lead 
to  nuaerloal  soheae  which  can  oonvargence  as  fast  as  the  current  error 
method. 

(5)  Along  the  line  of  designing  new  coaputatlonal  aathod.  a  systematic 
exploitation  of  tha  applicability  of  the  error  nathod  to  problaas  with 
elliptical  type  partial  differential  equation  nay  not  be  Inappropriate. 


APFBMSIX  k 


DAY 

Tuesday^  November  13»  1979  16:31 :02 


eiY  SPDl 

.FOR 

00100 

C 

SPDl— Sni.  UTION  FOR  LINEAR  HEAT  TRANSFER  PROBLEM. 

00200 

C 

ONh  DIMENRIONA!  UNSTFADY  HEAT  CONDUCTION 

00300 

n 

SUCCESSTVh  METHOD  SEPTEMBER  17.1979 

00400 

DIMENSION  U(lOl)  .UNdOl  )  .RSdOl )  .UNOdOl) 

0041)0 

OPEN  <  UN  I  T= 1 .  ACCESS^ '  SEQ I N '  .  D I  Al  OOUr.. ) 

00500 

READd.lO)  N.MAX 

00600 

10 

FORMAT (41 10) 

00700 

READd.20)  Ud).UI 

00800 

20 

FORMA T(RH 10.5) 

00850 

DY=1.0/FL0AT(N-1) 

00900 

DT=nY«»Y/2.5 

01000 

DO  30  1=2. N 

01100 

U(I)-=  HI 

01200 

30 

CONTINUE 

01300 

M=1 

01350 

TIMF==0,0 

01400 

WRITE<5.150)  N.NI  . T IMh . ( U ( I ) . 1  =  1 . N ) 

01500 

DO  40  1=1. N 

01600 

UN(I)=U(I) 

01700 

UNO(.T)=  UNd) 

01800 

40 

CONTINUE 

01900 

N1=N-1 

02000 

CFl=2.0)|;d  .0+DY«DY/DT) 

02100 

CF2=-  1 .0+DT/(DY#DY) 

02200 

50 

NL=0 

02300 

DO  60  1=2. N1 

02400 

IP=I+1 

02500 

11=1-1 

02600 

RS  <  I  >  :=  ( U  (IP ) +U  ( 1 1 )  -2 . 0«U  d  > )  /CF 1 +U  ( I  >  /CF2 

02700 

60 

CONTINUE 

02800 

70 

NT=0 

02900 

DO  90  1=2. N1 

03000 

IP=J+1 

03100 

11=1-1 

03200 

UN <  I  > :  (UN dP )  +UN (ID) /CF  1 +RR ( I ) 

03300 

DNO= ABS(UN(I)-UNO(I)) 

03400 

IF (DNO-0. 000001)  85.85.00 

03500 

80 

NT=1 

03600 

85 

UNO(  I  )-:UN(  I) 

03700 

90 

CONTINUE 

03800 

IF(NT-l)  120.100.100 

03900 

100 

NL=NI  +1 

04000 

00  TO  70 

04100 

120 

N=M+1 

04150 

TIME=TIMF+DT 

04200 

DO  130  1=1. N 

04300 

U(I)=UN(I) 

04400 

UNO(I).:UN(I) 

04500 

130 

CONTINUE 

04600 

WRITE(5.15())  M.NL.TIMF  .  (UN(  I )  .  1  =  1  .N) 

04700 

150 

FORMAT ( 1 OX . 21 1 0 . F 1 0 , 5 . /20 ( 5X . 6 ( F9 . 5 . 1 X ) / ) ) 

04800 

If (M-MAX)  50.200.200 

04900 

200 

STOP 

05000 

END 

APFBKDIX  B 


DAY 

Wednes:.(ji5</F  November  14 »  1979  li:37:07 


eiY  spii 

2.  FOR 

00100 

C 

SPD2— NFU  METHOD  FOR  LINEAR  HEAT  TRANSFERPROBLEMf 

00200 

C 

OHF  DIMENSIONAI.  UNSTEADY  HEAT  CONDUCTION 

00300 

C 

ERROR  METHOD 

00400 

DIMENSION  UCIOI.)  rUNdOl)  rRSdOl  )  rAdOl  )  »D<1.01. ) » 

00500 

IRdOl  )friFL(101) 

00600 

0PEN(UNIT-=1  tACCFSS-  'SEQIN'  »DJAL00UE) 

00700 

READddO)  NfMAX 

00800 

10 

FORMAT (4110) 

00900 

READ  d»20)  Ud>»UI 

01000 

20 

FORMAT  (8F).0*  5) 

01100 

DY=1.0/KL0AT(N-1) 

01200 

DT=:riY#DY/2.5 

01300 

DO  30  I=2»N 

01400 

U(  I  ):  =  UI 

01500 

30 

CONTINUE 

01600 

M=1 

01700 

TIME:=0.0 

01800 

WRITE(5rl50)  MrNI  dlME  »  ( U  d  )  f  I=d  r  N) 

01900 

DO  40  1=1 »N 

02000 

UN(I)=U(I) 

02100 

40 

CONTINUE 

02200 

N1=N-1 

02300 

CFl =2 . 0« (1 . 0+DY*DY/DT ) 

02400 

CF2=1 .0+DT/(DY»DY) 

02500 

50 

NL=0 

02600 

DO  60  1=2 »N1 

02700 

IP=I+1 

02800 

11=1-1 

02900 

RS  ( I )  =  ( U  <  I P )  +U  ( 1 1 )  -2  <  0)J:U  d  )  )  /CF 1 +U  (1 )  /CF2 

03000 

60 

CONTINUE 

03100 

70 

NT  =  0 

03200 

A(2)= 10,0**25 

03300 

E(2)=0.0 

03400 

DEL(N)-=0.0 

03500 

DO  130  I=2fNl 

03600 

IP=I+1 

03700 

11=1-1 

03800 

UNI=(UH(IP)+UN(I1 ) )/CHl+RS(I) 

03900 

DNO=ABS(UNI-UN(T)) 

04000 

Rd)=:UNT-UN(I) 

04100 

A(I  +  1  ):=CF1-)  .0/A(I) 

04200 

90 

B(I+1)=D(T)/A<I)-CF1*R(I) 

04300 

93 

IFdiNO. LE.0. 00001  )  GO  TO  130 

04400 

NT=1 

04500 

130 

CONTINUE 

04600 

IF(Nr-l)200fl40Fl40 

04700 

140 

DO  180  I=2fNl 

04800 

IR1=N-I+1 

04900 

IR2=IR1+1 

05000 

DEL ( INI ) ” DEL ( I R2 ) /A ( I R2 ) -B ( I R2 ) /A ( I R2 ) 

05100 

UNdRl  )==UN(IR1  )+DFJ.  (JRl) 

05200 

180 

CONTINUE 

05300 

NL=NI.  +  1 

05400 

GO  TO  70 

05500 

o 

o 

M=M+1 

05600 

TIME^dlMF+DT 

05700 

DO  230  I=1»N 

05800 

Ud)  =UNd) 

05900 

230 

CONTINUE 

06000 

WRITE(5»150)  MfNLfTIMEf  (UNd)»I=  1  »N) 

APFBNDIX  C 


DAY 

Tuesday »  Novt  uthct-  13 »  1979  19:09:36 
(»TY  SPD3.F0R 

00100  C  SPD3— DIFFUSION  IN  DOUBLE  TUBE  FLOW  PROBLEM f JULY/31/79 

00200  n  SUCCESSJVh  METHOD 

00300  DlMFK'SHm  T  ( 2  f  .1  0.1. )  r  YH  ( 2 » 101)  >  YO  ( 2 » 101 )  f  U  ( 2 , 101 )  ,  RO  ( 2 » 101 )  f 

00400  1C  ( 101 )  rDl.TdOi  )  fDI.TH(  101 )  rRO  I  (101)  »  TN(  101 )  »YFN(101 ) » 

00500  2Y0N(101 )  »ROYF(  101 )  ,ROYO(  101 )  f  DLYF(  101 )  »DYFN(  101 )  >- 

00600  3R(10J  )  fDl.YOdOl)  rDYONdOl) 

00700  OPEN  <  UNI T=:  1 1  AC(;ESS=  '  SEQI N  '  » DIALOGUE ) 

00800  READdrlO)  M » N » RF »  PR » SC  >  TW » WP » C  J 

00900  10  F0RMAT(2I10r lOH 10.5) 

01300  DR=0.5/FL0AT(N-1) 

01350  DX=DR»DR*RF/2.5 

01360  UR  I TECU  1  li )  M » N » DX  f  RE » PR  fSCfTUr  WP  f  C  J 

01370  15  F0RMAT(2I5»R(F10<5rlX)/) 

01380  C  BOUNDARY  CONDITIONS 

01400  Rd)=0.0 

01500  N1=N-1 

01600  DO  16  J--=1,N1 

01700  J1=J+1 

01800  R(  J1  )==R(  J)+DR 

01900  16  CONTINUE 

02000  LNI-0.5!eWP/DR+l  .0 

02100  C0P=8 . 0«DR)f:*2*RE*PR/DX 

02200  COS=CnF-*SC/F-R 

02300  DO  30  J=1,N 

02400  AJ=J-1 

02500  DO  25  1=1,2 

02600  T  ( I ,  J )  ^:  1 . 0+  <  I  W- 1 . 0 )  ft  ( A  J«riRft2 . 0 )  **2 

02700  R0(I,J)=l,0/T(Ir J) 

02800  IF(J-LNI)  17d7d9 

02900  17  YF(I,J)=  1.0 

03000  Y0(I,J)=-0.0 

03100  U( I , J)=CJftd .0-(2.0*R( J)/WP)**2) 

03200  GO  TO  25 

03300  19  YF(TrJ):  0,0 

03400  Y0(I»J)'0.5 

03500  U  <  I  f  J )  -•  1 . 0-  ( 2 .  OftR  ( J ) )  ft*2+  ( 1 , 0“WP**2 )  ftALOG  d  .  0/  ( 2 . 0*R  ( J ) 

03600  1)  )/AI.OG(Wr’) 

03700  25  CONTINUE 

03800  C  (  J )  ==  U  d  r  J )  ftRO  d  ,  J ) 

03900  30  CONTINUE 

04000  C(N)==0.0 

04100  C(2)=Cd) 

04200  C(LNI):-0.25*(C<LNI  +  1)  IC(I.NI-l)  ) 

04300  Ud,N)==0,0 


04400 

IJ  ( 2 » N )  ==  0 . 0 

04500 

Y0d,LNI)=-0.25 

04600 

Y0(2»LNI)-0.25 

04700 

YFdfLNI)=0.5 

04800 

YF(2fLNI)=  0.5 

04900 

IK=1 

05000 

SMTUR=0.0 

05100 

SMUR:-0.0 

05200 

DO  33  J=2,N1 

05300 

AJ1---J-1 

05400 

SMTUR  =  SMTUR  +  T ( 1 , J ) ftO (1 r  J ) ftA J1 

05500 

SMUR-SMUR+Ud  ,  J)«AJ1 

05600 

33  CONTINUE 

05700 

TBBB-=HMfUR/SMUR 

05800 

QW=260.0ft(TW-Td,Nl)  )/l.iR 

05900 


ANULT»QW/ <  260 , Oft ( f W- f BBB ) ) 


06000 

06100 

06200 

06300 

06400 

06500 

06600 

06700 

06800 

06900 

07000 

07100 

07200 

07300 

07400 

07500 

07600 

07700 

07800 

07900 

08000 

08100 

08200 

08300 

08400 

08500 

08600 

08700 

08800 

08900 

09000 

09100 

09200 

09300 

09400 

09500 

09600 

09700 

09800 

09900 

10000 

10100 

10200 

10300 

10400 

10500 

10600 

10700 

10800 

10900 

11000 

11100 

11200 

11300 

11400 

11500 

11600 

11700 

11800 

11900 

12000 

12100 

12200 

12300 

12400 

12500 


LINECCFK'fftAI,  Ff.C*) 


¥RiTF(J5»53;>)  IKfNITf  (T(l»  J)  ,  J=1  rN) 

WRITE (5» 570)  ANUITrOU 
WRITE(5»950)  < U ( 1 » J ) » J=1 » N ) 

WRITE ( 5 »  565 )  X »  < RO < 1 >  J ) , J= 1 , N ) 

WRITE ( 5  r  566  >  ( YE ( 1 , J ) f  J= 1 , N ) 

WRITE(l.r568)  (Y0(1»J)»J=1>N) 

X=0,0 

43  IK=IK+1 

IE(IK-?4)46f46rl703 
46  X=X+DX 

:  UNKK'OWS  NEAR  THE  CENTRAI  LINE(CFN  fRAI,  Ff.C*) 

NIT  =  1 

48  DO  57  1=1 f 2 
T(I,2)=T(Ifl) 

YF(If2)--YE<I»l) 

Y0(I.2)=Y0(I,1) 

57  CONTINUE 

;  UNKNOWNS  IN  THE  FIELD<2rN1) 

61  NT=0 

DO  1300  J=2»N1 

J1=J-1 

AJ1=J1 

AJ12^  <AJl+0.5)/AJl 
AJ32=(AJl-0«;>)/AJl 
CCOP==  COR«f;  ( J )  +4  ,  OJjtT  ( 1  r  J ) 

SQT2=SQRT  (  ( 0 . 25«C0P«(;  ( J )  > **2  }  COF*C  ( J )  *T  ( 1 ,  J) +0 , 5*A J1 2* 
1(T(1  r  J+l)+l  (:>r  J  +  1)  )!(:«:.>+0.5«AJ32>K(T(1»  J-1)+T(2»J~1)  )*)|c2) 
TN  ( J ) ::  -0 . 25«CC0P+HRT2 

TUS=T ( 1 » J+ 1 ) +T  <  2 » J+ 1 ) iT (1»J)+T(2»J) 

TLS=  T ( 1 r  J- l)+T(2»J-l)+T(lrJ)+T(2»J) 

COZJr-COS^C  (J)  I A  J1 2*TUS<  A  J323KTLS 

YFN(  J)  =  (AJ12>|:  iUS«<  YH  ( 1 ,  J+1 ) +YE  ( 2 »  J  +  1 )  ) +AJ32*TLS)t; 

1  (YF(1,  J-1  )  +  Vl  (.?fJ-l)  )  +  (C0S*C(J)-AJ12*TUS-AJ32>KTL.S)* 

2  YEO  »J))/(;02J 

Y0N(  J)::  (AJ12<:TUS!t!(  YOd  » J+1  )+Y0<2»  J+1 )  )  »AJ32)t!TLS* 

1  ( YO  ( 1  r  J- 1 )  +  YO ( 2  r  J- 1  )  )  +  ( ( J )  -A J 1 2*TUS~A J32*TLS )  X: 

2  YOdf  J))/C02J 
ROT( J)=TN( J)-l (2r J) 

ROYF(  J)=  YFN(  J)-YF  (2r  J) 

ROYO  <  J )  YON  ( J )  -  YO  ( 2 » J ) 

AROT=  ADS(ROT( J) ) 

AROYE"AHS(ROYh ( J) ) 

AROYO=ABS(ROYO(J)) 

IF<AR0T-0, 000001 )  572f572dl50 

572  IE ( AROYH -0.000001 >590 r 590 d 150 
590  IE ( AROYO-0 . 000001 ) 1300 » 1 300 » 1 1 50 
1150  NT=1 
1300  CONTINUE 

1F<NT)  1567d.567rl330 
1330  NIT=:Nn  +  l 

no  1538  J=2fNl 
T(2»  J)=^TN(  J) 

YF  <2fJ)-YhN<J) 

YO<?r  J)-=YnN(  J) 

1538  CONTINUE 
T(2fN)=TW 
YE(2»N)=YF(2f N1 ) 

Y0(2»N)=  Y0<2»N1) 

T(2d)=T(2r2) 

YE  ( 2  r  1 )  ==  YH  <  2  r  2 ) 

Y0<2f  1  )‘=Yn<2»2) 
no  TO  61 

1567  QW=260,0«:(U.'-T(2»N1)  )/i:(R 
DO  919  J=1»N 
R0(2f J)=  1.0/T(2f J) 

U<2» J)=  n( J)/R0(2» J) 


1567 


26 
12700 
12800 
12900 
13000 
13100 
13200 
13300 
13400 
13500 
13600 
13700 
13800 
13900 
14000 
14100 
14200 
14300 
14400 
14500 
14600 
14700 
14800 
14900 
15000 
15100 
15200 
15300 
15400 
15500 
15600 
15700 
15800 

e 


CONTINUE 
SMTUR-O.O 
SMUR:  0.0 
00  908  J=2»N1 
J1=J-1 

SMTUR=SKTIIR+  f  <  2 1 J )  «0  ( 2  r  J )  #R  ( J )» ( 2 , 0*DR ) 
SMUR=SMUR+U ( 2  f  J ) «R ( J ) « ( 2 . 0*DR ) 

CONTINUE 

TBBB^SMTUR/SM'JR 

ANULT=QW/  ( 260 . 0^<  ( TW-TBBB )  ) 

WRITE 535)  IKf NITf (T(2f J>  f J=1^N> 

WR I TE ( 5  r  570  >  ANUL T  r  QW 
WRITF.<5»950)  ( U ( 2 » J )  » J=1 » N ) 

WRITE<5»565>  X » ( RO ( 2 » J ) » J=1 » N > 

WRITH (5r566)  ( YF(2f J) f J=1 yN) 
WRITE(5y568) (Y0<2»J)»J=1>N) 

F0RMAT<2Xr '  I yOW' y2I5»5(F9.5f lX)/25 
I  (17X»5(F9.5ylX)/)) 

FORMAT ( 1 7X  r  F9 . 5 » 1 X  y  F9 . 3 ) 


1679 
)  697 
1703 


FORMAT  (2Xr‘'  U'yl 
F0RMAT(2Xy'  XyRO'.h 
F0RMAI(2Xr'  YF ' y J 
F0RMAK2Xr'  YO'yl 
»()  1679  J=lyN 
T<ly  J):-T(2y  J) 
YF(lyJ)==YF(2yJ) 

Y0(1 yJ)-  Y0(2y J) 

R0(]  yJ)::1.0/T(lyJ) 

U<lyJ)=C(J)/Rn(ly J) 

CONTINUE 

GO  TO  43 

STOP 

END 


U'  yl0Xy5(F9.5ylX)y/25(17Xy5(F9.5ylX)/)  ) 

X  y  RO ' »  H  0 . 6  y  5 ( F9 . 5  y 1 X ) y /25 ( 1 7X  y  5  <  F9 . 5  y 1 X ) / ) ) 
YF' y  J0Xy5(F9.5ylX)/25(17Xy5(F9.5y  lX)/) ) 

YO' yl0Xy5(F9,nylX)/25<17Xy5(F9.5ylX)/) ) 


APPENDIX  D 


DAY 

T  uesdau 
eiY  SPD4 
00100 
00200 
00300 
00400 
00500 
00600 
00700 
00800 
00900 
01300 
01350 
01360 
01370 
013R0 
01400 
01500 
01600 
01700 
01800 
01900 
02000 
02100 
02200 
02300 
02400 
02500 
02600 
02700 
02800 
02900 
03000 
03100 
03200 
03300 
03400 
03500 
03600 
03700 
03800 
03900 
04000 
04100 
04200 
04300 
04400 
04500 
04600 
04700 
04800 
04900 
05000 
05100 
05200 
05300 
05400 
05500 
05600 
05700 
05800 
05900 


f  NovcinlxM' 
.FOR 


1979  i9:o8:oo 


3Pri4--»IFHJft.T0R  JN  DOUBLE  TUBE  FLOW  PROBLEMf JULY/31/79 
ERROR'  METHOD 

DIMENSION  T<2»101 >  rYF(2f 101 ) »Y0(2»101) fU(2f 101) fR0(2f 101 ) f 
1 C  ( 1 01 )  f  1  N  ( 3  01  )  r  Y)-  N  ( 101 )  f 

2Y0N < 1 01 >  r  A7  < 1 01 ) » BT ( 1 01 ) f  A Y ( 1 01 ) f  BYF ( 1 01 ) » 

3R(101)  »BY0<10i)  r);i:;i.  (101)  »DEF(101)»DE0(101) 

OPEN  ( UNIT':  1  r  AOOi; SS^  ' SEQIN '  x  DIALOGUE ) 

READ (1x10)  M  X  N  X  Rr  X  PR  X  SC  X  TW » WP  x  C J 
3  0  FORMA  / ( 2 X 1 0  x 1 OF 1 0 . 5 ) 

DR=0.ri/H.OAC(N-l) 

DX=DR*DR*RF/2.5 

WRI TF ( 5  X 1 5 )  M  X  N  X  DX  X  RE  X  PR  X  SC  X  TW  X  WP  X  C J 

15  FORMA  I  ( 2;i;r) x  O  <  no .  5 x  l  X )  / ) 

BOUNDARY  CONDITIONS 
R(1)=0.0 

N1=N--1 

DO  16  J=lxNl 
J1=J+1 

R(  Jl)==-R(  J)  +  DR 

16  CONTINUE 
LNI=0.5*UP/DR+1 .0 
C0P=8 . 0*DR«<:2«Rh  «PR/DX 
COS^  COP#S(;/PR 

DO  30  J=lxN 
AJ=J-1 
DO  25  1=1x2 

T  ( I X  J )  =  1 . 0+  ( TW- 1 . 0 ) « ( A  J«DR^:2 . 0 )  J(:«2 
R0(IxxJ)==1.0/T(l,j) 

IF(J-I..NI)  3  7x17x19 

17  YF(IxJ)=1.0 
Y0( I X J)”0.0 

U  ( I X  J  )  ••=  C  J#  ( 1 . 0-  ( 2 . 0*R  ( J  )  /WP )  *«2 ) 

00  TO  25 
19  YF(IxJ)=0.0 
Y0(IxJ):=0,5 

U ( I  X  J ) :  1 , 0- <  2 . 0*R ( J ) ) #«2+ ( 1 . 0-WP**2 ) *ALOG ( 1 . 0/ ( 2 , 0*R ( J ) 

1) )/ALOG(WP) 

25  CONTINUE 

C  <  J )  =•  U  ( 1 X  J )  sKRO  ( 1 X  J ) 

30  CONTINUE 
C(N):  0.0 
C(2)=  (;<  1 ) 

C(LNI)-  0.25# (0(1. NI  +  1  )TC(LNI-1)  ) 

U(lxN)=  0.0 

U(2xN)=0.0 

Y0(3  xl.N.1  )>=(). 25 

Y0(2xl.NT):=0.25 

YFdxl.Nl)  0.5 

YF(2xl..Nl)=:0.5 

IK=1 

SMTUR=0.0 
SMUR=0.0 
DO  33  J=2xNl 
AJ1=-J-1 

SMTUR=  SM  f  UR+  T ( 1 x  J ) #U ( 3  x  J ) #A J1 
SMUR=SMUR+U ( 1 X  J ) #A J1 
33  CONTINUE 

TBBB-SMTUR/SMUR 

QW=260. 0# ( TW-T ( 1 x  N1 ) )/DR 

ANUl.T"  «W/  ( 260 , 0#  (  I W-TBBB )  ) 


06100 

06200 

06300 

06400 

06500 

06600 

06700 

06800 

06900 

07000 

07100 

07200 

07600 

07400 

07500 

07600 

07700 

07800 

07900 

08000 

08100 

08200 

08300 

08400 

08500 

08600 

08700 

08800 

08900 

09000 

09100 

09200 

09300 

09400 

09500 

09600 

09700 

09800 

09900 

10000 

10100 

10200 

10300 

10400 

10500 

10600 

10700 

10800 

10900 

11000 

11100 

11200 

11300 

11400 

11500 

11600 

11700 

11800 

11900 

12000 

12100 

12200 

12300 

12400 

12500 


WRITE<5»53S)  IK»NTT» (T(1/J)»J=1»N)  ' 

WRITF(t5»;)70>  ANULT»QW 
WRITF<'ir950)  < U<  1  f  J) » J=1  f  N) 

WRITE (5» 565)  X» (R0( 1 » J) » J-1 ,N) 

WRITE(;m5"66)  (YF(1»  J)  »  J=1»N) 

WRITE(5»56B)  ( Y0< 1 » J) » J=1 ,N> 

X=0.0 
43  IK=IK+1 

IF(IK-M)46r46»1703 
46  X=X+DX 

C  UNKNOWS  NEAR  THE  CENTRAI  l.INF((;HN  (RA!  B.C.) 

NIT=1 

48  no  57  1=1 » 2 
T  ( I » 2 )  =•  I  <  I » 1 ) 

YF  ( .T  r  2 )  ==  YK  <  I » 1 ) 

YO  ( I  f  2 )  =  Yo  ( ;i;  r  1 ) 

57  CONTINUE 

(;  UNKNOWNS  IN  THE  FIELD(2»N1) 

61  Nf=0 

AT (2) =  1.0 
BT(2)==0.0 
WO  1300  J=2fNl 
J1=J--1 
AJl^Jl 

AJ12^(AJl+0.5)/Ajl 
AJ32-  (AJl-0.:))/AJl 
CCOP=- COPCC  ( J )  +4 . 0*T  ( 1  r  0 ) 

SQT2=SQRT  ( <  0 . 25>»!C0P*C  <  J ) )  **2+C0P*C  (J)*T(lfJ)+0. 5*A  J1 2 
1  ( T  ( 1  r  JH  )  +  I  ( 2  r  J  {•  1 )  )  «*2+0 . 5* A  J32»  ( T  ( 1 » J  - 1 )  I  T  ( 2 » J  - 1 )  )  *» 

TN  ( J )  ==  -0 . 25*CrOP+SQT2 
R0T=TN<J)-T(2» J) 

AROT=ABS(ROT) 

CAT=2.0)»:SQT2/(AJ125f:(  (  (1,JH)+T(2»J+1)  )  ) 

CBT^  AJ62Je:(  f  (Ir  J-1  >  +  1(2,  J-1 )  )/( AJ12*(T(  1 » J+1  )  +  T<2f  J+1 )  )  ) 
AT(  J+1  )=•  CAT-rBT/A/f'J) 

BT  <  J+ 1 ) :  :CWr«BT ( J ) /AT ( J ) -CAT*ROT 
345  IF(AR0T~0, 000001 )  1 300 , 1300 » 1 150 
1150  NT=1 
1300  CONTINUE 

IF(NT,LE.O)  GO  TO  1567 
DEL(N3 )  -WT(N)/Af  <N) 

T(2fNl)=  •f(2»Nl)TWEL(Nl) 

WO  3313  J=2fNl 

NJ=N-J 

NJ1=  N+l-J 

nEL(NJ)-:  (Wh.L(NJl)-BT(NJl)  )/AT<NJl) 

T<2»NJ)-  T(2»NJ)  +  WI-:i.  (NJ) 

1313  CONTINUE 

T(2»l)=-T(2»2) 

T(2»N)=TW 
NIT=N7T+1 
00  TO  61 
1567  NY=0 

AY(2)=1.0 
BYF(2):^  0.0 
BY0<2)=0.0 
WO  1400  J=2»N1 
J1=J-1 
AJ1=J1 

AJ12-  (AJl+0,5)/AJl 
AJ32- (AJl-0,5)/AJl 

TUS=T<1» J+l)+T(2» J+1)+T<1» J)+T(2» J) 

1LS=  T<l»U)+f  <2»  J)n  (If  J-1  )+T(2»  J-1) 

C02 J  (;0SK:0  ( U )  +A J12*TUS+A J32*TLS 

YFN(  J):  (AJ3  2#  fOH)i<(YF(lf  J+l)+YF(2f  J+1)  )+AJ32*TLS# 

1  JYFdr.J  l)  +  YK2r.J-i))  +  <C0S#C<J>-AJ12*TUS-AJ32#TLS)* 


P 


ro  * 


12600 
12700 
12800 
12900 
13000 
13100 
13200 
13300 
13400 
13500 
13600 
13700 
13800 
13900 
14000 
14100 
14200 
14300 
14400 
14500 
14600 
14700 
14800 
14900 
15000 
15100 
15200 
15300 
15400 
15500 
15600 
15700 
15800 
15900 
16000 
16100 
16200 
16300 
16400 
16500 
16600 
16700 
16800 
16900 
17000 
17100 
17200 
17300 
17400 
17500 
17600 
17700 
17800 
17900 
18000 
18100 
18200 
18300 
18400 
18500 
18600 
18700 
18800 
18900 
19000 
LOO 


139: 


1400 


YON (  J ) :  ( A J.1 2*TUS*  ( YO <  1 ,  J+ 1 )  +  YO <  2 » J+ 1  >  )  +A  J32*TLS* 

(Y0<1»  J-1)+Y0(2f  v)-l)  )  >  (r.OS)»cC(  J)-AJ12*TUS-AJ32!tcrLS>* 
Y0<1»  J)  )/ClJ2J 
ROyF=YFN<J)-Yf-  <2fJ) 

R0Y0=-Y0W<J)-Yn(2»J) 

AROYF==AHS(ROYF) 

AROYO^  AEiS(ROYO) 

CAY=nt)2J/(AJ12*TUS) 

CBY=A  J32)f:TLS/  ( A  J1 2*TUS ) 

AY  ( J+1 )  ==  f;AY-f;BY/AY  ( J  ) 

BYF  <  J+ 1 ) =CBY*B YF ( J ) /AY ( J ) -CAY*ROYF 
BYO  ( J+1 ) :  IIBYCBYO ( J )  /AY  ( J )  -CAY»ROYO 
IF (AROYF.LF.O* 000001)  GO  TO  1395 
K'Y=1 

IF (ARCiYO.  I  F,  0.000001)  GO  TO  1400 
NY=1 

CONTINUE 
IF(NYJ.F.O)  GO  TO  1576 
DEF ( N1 ) =BYF ( N ) / ( 1 . 0-AY ( N ) ) 

DE0(N1 ): BYn(N)/(1.0-AY(N) ) 

YF(2»N1 ):  YF(2»N1)+DEF(N1 ) 

YO ( 2 » N1 ) - Y0<  N1 ) +DEO  <  N1 ) 

):iO  1414  J=2»N1 

NJ=N-J 

NJ1=N+1-J 

DEF(NJ):=  0:iFF<NJl)~BYF(NJl)  )/AY(NJl) 

DEO ( N J )  ( BE  0  ( N J1 )  -BYO ( N J 1 )  )  /AY  ( N J 1 ) 
YF(2»NJ)=YF<2fHJ)+IiEF(NJ) 

YO  <  2 » N  J  )  YO  (  2 » N  J  )  +DEO  ( NJ ) 

1414  CONTINUE 

YF(2fl)==  YH  (.?»?) 

YF(2fN)-  YE  (,?rNl> 

Y0(2rl  ):-Y0(2f2) 

Y0(2»N)=Y0(2»N1) 

NIT=NIT+1 
GO  TO  1567 

1576  nU=260.0«(TW-T(2»Nl) )/DR 
I.IO  919  J=lfN 
R0(2»  J)-^1.0/T(2»  J) 

U(2»  J)=-C(  J)/R0(2»  J) 

919  CONIINUE 
SMTUR=0.0 
SMUR=0,0 
)){)  908  J=2fNl 
J1=J-1 

SMTUR^  BM  rUR+T  ( 2  r  J )  «U  <  2 ,  J )  )i:R  ( J ) « <  2 . 0*DR ) 

SMUR  SMUR+U  ( 2  r  J )  «R  <  J ) )« ( 2  ♦  0*DR ) 

908  CONTINUE 

TBBB  8K(UR/KMl.lR 

ANULT==  nw/ <  260 , 0)f:  ( TW-TBBB )  ) 

WRITE(5»535)  IK»NIT»  <T(2f J) f J=1 »N) 

ANULTfOW 
(U(2»J)»J=1,N) 

X»(R0(2»J)»J=1»N) 

(YF(2f J)» J=1»N) 

WRITE (5r 568) <Yn<2»J) »J=1»N) 

535  FORMA l(2Xf'  Tf OW ' » 2I5f 5(F9. 5f lX)/25 
1  (17X»5<F9.5»1X)/)) 

570  FORMA! <17X»F9.5»lXfF9.3) 

950  F0RMAT(2X»'  U ' » 1  OX » 5 (F9. 5f IX ) » /25 ( 1 7X» 5 <F9 , 5» IX )/ ) ) 

565  F0RMAT(2X»  '  X»RO'  rHO.  6r5(F9,5f  IX)  »/25(  17X»5<F9,5»  IX)/)  ) 

566  F0RMAf(2X»'  YF  '  ,  lOX » r.(F9 . 5 » IX ) /25(  17X » 5(  F9 .5»  IX  )/>  ) 

568  F0RMAT(2Xf'  YO ' » 1  OX r 5( F9 , 5» IX ) /25 < 17X »5(F9. 5 dX )/ ) ) 

BO  1679  J»1»N 
TUyJ)«T<2r  J) 


WRITE (5f570) 
WRITF<r.r950) 
WRITE (5» 565) 
WRITE(5»566) 


19200 

YF<1,J)»VF<2»J> 

19300 

Y0<1»  J)=-  Y0(2»  .1) 

19400 

R0(1»J)=  1.0/f (1» J) 

19500 

U(l»  J)==(MJ)/K()(lr  J) 

19600 

1679 

CONTINUE 

19700 

1697 

GO  TO  43 

19800 

1703 

STOP 

19900 

END 

e 


APPENDIX  E 


DAY 

TuesdiMif  November  13  r  1979  201 03  J 13 
etY  SPD5.F0R 


00100 

C 

SPD5— NON  LINEAR  TUBE  FLOW  HEAT  TRANSFER  PROBLEM f 

00150 

C 

SUCCESSIVE  METHOD  NOVEMBER  1»1979 

00200 

DIMENSION  T(2»101  ) »  U( 2»  101 ) , RfJ ( 2 r  101 )  » 

00300 

ICdODrTNdODt 

00500 

3R(101) 

00600 

OPEN  ( UN  J T=:  1  f  Af;nFSS==  ' SF.  R IN ' » B I  Al  OSJ IE > 

00700 

READ<1»10)  M»N»R!- »PRrTW 

00800 

10 

FORMAT (211  Or  lOHO. 5) 

01200 

DR=0.;)/Hl.nAT<N-l) 

01250 

DX=DR«riR)f;RE/2.5 

01260 

WR I TE ( 5  r 1 5 )  M  »  N  r  DX  r  RF  r  PR  r  T W 

01270 

15 

F0RMAT(2JrM8(F10.5f  IX)/) 

01280 

C 

BOUNDARY  CONDITIONS 

01300 

R(l)=  0.0 

01400 

N1=N-1 

01500 

DO  16  J=lrNl 

01600 

J1=J+1 

01700 

R( J1 )=R( J)+DR 

01800 

16 

CONTINUE 

02000 

C0P=:8 . 0«DR«»23»tRE»PR/DX 

02200 

))D  30  J=lrN 

02300 

AJ==J-1 

02400 

))0  25  1  =  1  r2 

02500 

T  ( I  r  J )  =^1 . 0+ ( TW- 1 . 0 )  *  ( A  J*DR*2 . 0 )  !Sc!f:2 

02600 

R0(I»J)=1.0/T(J»J) 

03000 

Udr  J):  1.0-(2.0«R(J))*>5;2 

03600 

25 

CONTINUE 

03700 

C( J)  ■^Udr  J)*RO(lr J) 

03800 

30 

CONTINUE 

03900 

C  ( N )  ••=  0 . 0 

04000 

■UdrN)=0.0 

04100 

U(2»N)=0.0 

05200 

IK=1 

05300 

SMTUR==^0.0 

05400 

SMUR=:-0.0 

05500 

DO  33  J=2»N1 

05600 

AJ1=J-1 

05700 

SMTUR=SMTUR+T  d  r  J )  «U  d  r  J )  JScA  J1 

05800 

SMUR:=SMUR+Ud  r  J)«AJ1 

05900 

33 

CONTINUE 

06000 

TBBB^^^SMfUK/SMUR 

06100 

QW=260 ♦ 0# ( TW-T ( 1 r  N1 ) ) /DR 

06200 

ANULT==nW/  ( 260 .  OC  <  TW-TBB)i ) ) 

06300 

URITF(r;r535)  IK  r  NIT  r  (  Td  r  J  )  r  J-=  1  r  N ) 

06400 

WRITE(5»570)  ANULTrQU 

06500 

WRITE<5»950)  (Ud  r  J)  r  1  r N) 

0o600 

WRITE(r>r5A5)  Xr  (ROdr  J)  r  J=:l  rN) 

06900 

X=0.0 

07000 

43 

IK=IK+1 

07100 

IF(IK-M)46r46rl703 

07200 

46 

X=X+DX 

07300 

C 

UNKNOWS  NEAR  THE  CENTRAI.  LINE ( CEN 1 RA!.  B.C.) 

07400 

NIT=1 

07500 

48 

DO  57  I=lr2 

07600 

T(I»2)»T(I»1) 

07700 

U<I»2)-U(Irl) 

07800 

R0<Ir2)=R0(Irl) 

08100 

57 

CONTINUE 

08200 

C 

UNKNOWNS  IN  THE  FIELD(2rNl) 

08300 

61 

NT=0 

08400 

08500 

08600 

08700 

08800 

08900 

09000 

09100 

09200 

10100 

10400 

10700 

11000 

11100 

11200 

11300 

11400 

11700 

11800 

12250 

12300 

12400 

12420 

12440 

12460 

12480 

12500 

12600 

12700 

12800 

12900 

13000 

13100 

13200 

13300 

13400 

13500 

13600 

13700 

14000 

14100 

14200 

14300 

14400 

14700 

14800 

15100 

15200 

15300 

15400 

15500 

15600 


1 1 50 
1300 


1313 


1567 


1679 

1697 

1703 


JlO  1300  J*2»N1 

J1=J-1 

6J1=J1 

6J12=^<AJl+0.5)/AJl 
6J32=^  (AJl-0,5)/6Jl 
CCQP=^C0P»C(J)  +  4«0#T  (1,  J> 

SQT2=SQRT  ( <  0 . 25«f;fJP«(;  ( J  > )  «*2+C0P*C  <  J ) « T  ( 1  r  J )  +0 , 5*AJ12J|t 
.  <  T  ( 1  r  J+ 1 )  +1  ( 2  r  J+ 1 )  >  *»2+0 . 5«A.I32«  ( T  ( 1.  r  J  - 1 )  +  T  <  2  f  J- 1 ) )  3»c»2 ) 
TN( J)=  -0.25»CC0P+SQT2 
ROT=TN( J)-T(2» J) 
rtROT=^ABS<ROT) 

IF  <  AROT-0 « 000001 )  1 300 » 1 300 » 1 1 50 
NT=1 

CONTINUE 

IF<NT.I.E,0)  GO  TO  1567 
no  1313  J=2»N1 
T(2»  J):  TK'(  J) 

CONTINUE 
T(2»l  )=  !  (?»  2) 

NIT=NIT+1 
GO  TO  61 

GW==260.0«(TW“T(2»N1  )  >/)iR 
no  919  J=1»N 
R0(2»  J)-^  1  .0/1  (2»  J) 

U(2f  J):a;<  J)/K0(2r  J) 

CONTINUE 
SMTUR':0.0 
SMUR=0,0 
DO  908  J=2fNl 
J1=J-1 

SMTUR=SMTUR+  T  <  2 » J ) «U  <  2  r  J ) «R  <  J ) * ( 2 . 0*DR ) 

SMUR=SMUR+U  ( 2 » J )  «R  <  J ) « ( 2 . 0)S:):tR ) 

CONTINUE 

TBBB*SNTLIK/BMUR 

ANULT=QW/ ( 260 . 0« ( ( W-TBBB ) ) 

WRITF(!)f53o)  IK»NITf  (T(2f  J)  f  J  :1,N) 

URITF(5»570)  ANULT»QW 
WRlTF(li»950)  (U(2»  J)  f  1  rN) 

WRITE (5» 565)  Xf ( R0(2 » J> » J -  1 r N) 

F0RMAT(2X» '  T »QW ' » 2I5» 5 ( F 9 . 5 » IX ) /25 
I  (17Xf5<F9,5»lX>/)) 

FORMAT ( 1 7X  r  F9 . 5  F 1 X » F9 . 3 ) 

F0RMAT(2Xf '  U' f 10Xr5(F9.5F IX) »/25( 17XF5<F9.5r IX)/) ) 
FORMAT <?Xf  '  XfRO'  FF10.6F5(F9,5rlX)  r/2:;(l/Xf5(F9.5FlX)/)) 
»()  1679  J=1fN 
T(1f J)= T<2» J) 

R0<1f  J)==1.0/T(1f  J) 

U(1f  J)==C(J)/R0(1  F  J) 

CONTINUE 
GO  TO  43 
STOP 
END 


AFfBHSIX  F 


DAY 

Tuesday*  November  13*  1979  20:28:37 
eiY  SPD6.F0R 

00100  C  SPD6--N0NLINEAR  TUBE  FLOW  HEAT  TRANSFER  PROBLEM 

00150  C  ERROR  METHOD*  NOVEMBIK  4*1979 

00200  DIMENSION  T(2*  101 )  *11(2*  101 )  *R0(2*  101 )  * 

00300  1C(101)*A(101) *B(101)*TN(10i)* 

00500  3R(101>*DFL(J01)*Tn(2*10J ) 

00600  OPEN(HNIT:=l»Af;(;lr  SS-  '  SFQIK' '  *  DIAI  OOHE ) 

00700  READdrlO)  M*N*RFfPR*TW 

00800  10  FORMA) (2110* lOFlO. 5) 

01200  DR=0.  r//FI.  DA  r(N-l ) 

01250  DX=:DR#DR!ttRE/2.5 

01270  IP=1 

01275  WRITF(5*15)  M* N * DX * RF *  PR  *  1 W 

012R0  15  F0RMAT(2I5*8(F10.5*1X)/) 

01285  C  BOCiWDARY  CONDITIONS 

01300  R(l)-0.0 

01400  N1=N-1 

01500  DO  16  J==1*N1 

01600  J1=J+1 

01700  R(J1  )-R(>J)  +  riR 

01800  16  CONTINUE 


02000 

02200 

02300 

02400 

02500 

02600 

03000 

03600 

03700 

03800 

03900 

04000 

04100 

05200 

05300 

05400 

05500 

05600 

05700 

05800 

05900 

06000 

06100 

06200 

06300 

06400 

06500 

06600 

06900 

07000 

07050 

07100 

07200 

07300 

07400 

07500 

07600 

07700 

07800 

08100 


C0p=-  B .  0«);iR««2)!:RF«PK/DX  / 

):iD  30  J=1*N  / 

AJ=J-1 

DO  25  1=1*2 

T ( I  *  J ) ^  1 . 0+  (  TW-1 , 0 ) « ( A J*DR«2 ♦ 0 )  **2 
RO(I» J)=1.0/T(I*J) 

U(I*  J)  =  1.0-(2,0)5:R(  J)  )«#2 
25  CONTINUE 

C(  J)=Ud  *  J)«K0(1*  J) 

30  CONTINUE 
C(N)=0.0 
U(1*N)=0.0 
U(2*N)"0,0 
IK=1 

SMTUK  0*0 
SMUR=0.0 
DO  33  J=2*N1 
AJ1=J-1 

SMTUR- SM  f UR+T ( 1  *  J ) *U ( 1  *  J  >  «A J1 
SMUR=SMUR+U(1 * J)«AJ1 
33  CONTINUE 

TBBB=SMTUR/SMUR 

QW=260 , 0#  (  I  W-T  ( 1  *  N1 )  )  /))R 

ANULT  =  NW/ ( 260 ♦ 0# ( TW-TBBB ) ) 

WRITE  ( 5  *  535 )  IK  *  NIT  *  ( T  ( 1  *  J )  *  1 » N ) 

WRITE (5* 570)  ANULT *QW 
WRI TF  ( 5  *  950 )  ( U  ( 1  *  .J )  *  1  *  N ) 

WRITE ( 5  *  565 )  X  * ( RO ( 1  *  J ) *  J  1  *  N ) 

X=0.0 

43  IK=IK+1 

IP=IP+1 

IF(IK-M>46*46*1703 
46  X=X+DX 

C  UNKNOWS  NEAR  THE  CENTRA).  LINE(CENTRAL  B.C.) 
NIT=1 

DU  57  1=1*2 
T(I*2):  T(I*1) 

U(I*2)”U(1*1  ) 

R0(I*2)'=RfJ(I*l> 

57  CONTINUE 


08200 

08300 

083-10 

08380 

08400 

08500 

08600 

08700 

08800 

08900 

09000 

09100 

09200 

10100 

10400 

10420 

10440 

10460 

10480 

10700 

11000 

11100 

11200 

11220 

11240 

11300 

11350 

11400 

11450 

11500 

11700 

11800 

11820 

12250 

12300 

12400 

12410 

12420 

12430 

12440 

12500 

12600 

12700 

12800 

12900 

13000 

13100 

13200 

13300 

13320 

13340 

13400 

13500 

13600 

13700 

14000 

14100 

14200 

14300 

14400 

14700 

14800 

15100 

15200 

15300 

Afi4P0 


C  UNKNOUNS  IN  THE  TiEL  1)12 r  Nl  > 

61  NT=0 

A(2):=l,0 

B<2)*0.0 

WO  1300  J=2»N1 

J1=J-1 

6J1=J1 

AJ12"  (AJl+0♦t>)/A^)l 
AJ32=  (AJl-0i5)/AJl 
ccop  ^  ( j  > +4  *  o«T  ( 1 » j ) 

SQT2=SQRT  (  ( 0 . 2r>«(:()P#(:  ( J ) )  *«2+(:()P#f;  ( J )  «T  <  J  r  J )  +0 . 5»A  J12)t£ 

1  ( T  <  1  » J+  J  )  +  I  <  2  r  J  1 1  > )  #«2+0 . 5«A .132#  ( T  ( 1  f  J- 1 )  +  T  ( 2  r  J- 1 ) )  ##2 ) 
TN  ( J ) :  -0 . 25#(;f;()P+8nT2 
ROT=TN( J)-T(2f J) 

AR0T:=AB8(R0T) 

CA=2.0#SRT2/(AJ12*(T(1» Jfl)+T(2f J+l>  ) ) 

CB=A J32# (  f  ( 1  r  J- 1 )  +  T ( 2 » J- 1 ) ) / ( A J1 2* < T ( 1  r  JH  )  +T ( 2 » J+l ) )  ) 
A(J+1 )=  CA-CB/A( J) 

B  ( J+ 1  )  =-CP#B  ( J )  /A  ( J )  -CARROT 
IF ( AROT-O ♦ 000001 )  1300 » 1300 » 1 1 50 
1150  NT=1 
1300  CONTINUE 

IF(NT.LH<0>  GO  TO  1567 
DEL(?^1  )”-B<N)/A(N) 

T<2»Nl)=T(2»Nl)+riHI_<N1  ) 

no  1313  J=2»N1 

NJ=N-J 

NJ1=N+1-J 

DEL<K'J)=  (ItEI  (NJ1)-B<NJ1)  )/A(NJl) 

T  ( 2 » N  J )  =-T  <  2  f  N  J )  +IiH  I..  ( N  J ) 

1313  CONTINUE 

T(2»l):  T(2»2) 

T<2fN)-=TW 
NIT=NIT+1 
(50  TO  61 

1567  QW=260.0«(TW-T(2»N1) )/DR 
DO  919  J=1»N 
R0(2f  J)=-1.0/T(2r  J) 

U<2j  J)--f;(J)/R0(2r  J) 

919  CONTINUE 
SMTUR=  0.0 
SMUR=0.0 
DO  908  J=2fNl 
J1=J-1 

SMTUR=-SMTUR  +  T  ( 2 » J )  #U  ( 2  r  J )  #R  ( J )  #  (  2 , 0#DK ) 

SMUR- 8MUK+U (2»J)*R(J)*(2. 0#nR ) 

908  CONTINUE 

TBBB=SNTUR/SMUR 
ANULT=-(^U/(260.0#(TW-TBBB)  ) 

IF( IP.NH.il)  GO  TO  731 
IP=1 

URITF(5rr>35)  JK » NIT  »  (  (  ( 2 »  J ) » J=-l  f  N ) 

WRITE(5f570)  ANULTfOW 
WRITE(5»950)  (U<2» J) f J=1 *N) 

WRITE  ( 5 » 565 )  X » <  RO ( 2 r  J )  r  1  r  N ) 

535  FORMAT <2Xr'  T » QW ' » 2T5 » 5 ( F9 . 5 f 1 X ) 725 
1  (17X»5(F9.5»1X)/)) 

570  FORMAT ( 1 7X  f F 9 . 5  r IX » F  9 . 3 ) 

950  FORMAT ( 2X » '  U ' » 1 OX » 5 ( F9 . 5 » 1 X ) r  725 ( 1 7X  r ! )  ( F  9 . 5 » 1 X ) 7 ) ) 
565  FORMAT (2Xr '  X » RO ' » FIO . 6 r 5< F 9 < 5 » IX ) » 725( 1 7X » 5(F9. 5» IX ) 7 ) ) 
731  DO  1679  J=lfN 
T(l»  J)=-T(2»  J) 

R0(1»  J)=-i,07l  (If  J) 

U<l»J)==C(J)7K0<lr.l) 

1679  CONTINUE 


▲PFBNDIX  G 


DAY 

rije<:,d<)w t  Novembcv  13»  1979  21:1-1:26 


eiY  SPD7 

.F 

C^TY  SPD7 

.FOR 

00100 

C 

SPD7— HEAT  TRANSFER  UITH  RADIATION>AN  UNSTEADY? 

00200 

C 

ONF  DIKENS.TONAI  FIN  PROBI. FM r NOOh MKF R  13?  1979 

00300 

C 

SUCCESSIVE  METHOD 

00400 

DIMENSION  T(2?101)  rDdOl)  ?TN(101> 

00500 

OPEN(llNJT-=  1  r  ACCESS-:  ' SFQIN '» DIALOGUE ) 

00600 

READ  < 1 ? 1 0 )  GC  ?  CC  r  RC  r  ( R » B 

00700 

READ(1?20)  M?N 

00800 

10 

F0RMAT(8F10.5) 

00900 

20 

FORMAT (41 10) 

01000 

N1=N-1 

01100 

riX=:GC/FI  0AT(N1 ) 

01200 

DT=DX)f:riX/2.5 

01300 

WR I TE  ( 5  ?  30 )  GC  ?  CC  r  RC  r  TK  r  B  r  DX  ?  D  f  rl-UN 

01400 

30 

FORMAT ( 1 OX  ?  7F 1 0 . 5  f  2 J 1 0/ ) 

01500 

DO  50  I--::l?2 

01600 

DO  40  J=1?N 

01700 

T ( I »  J )  :  1.0 

01800 

40 

CONTINUE 

01900 

50 

CONTINUE 

02000 

IK  =  1 

02100 

NIT  =  1 

02150 

TIME=0.0 

02200 

WRITE(5»53)  IK?NJTrTIKiK  r  O  (Ir  J>?J=1?N) 

02300 

53 

F0RMAT(2X?2I5?F8.5? IX r 5 < F9 . 5 i IX ) / ? 20 ( 21X » 5 < F9 . 5 ? IX ) / > / ) 

02400 

56 

NIT=1 

02500 

58 

NT  =  0 

02600 

TNa)-T(lrl) 

02700 

DO  80  J-2?N 

02800 

TM=0.5«:(  1  (1  f  J)+T(2?  J)  ) 

02900 

D  ( J )  =  1  . 0+RC>!;  (  TM««3  +  4  .  ()«TR*TM**2+6 . 0*TR*)K2*TM+4 . 0)KTR)*(»3 ) 

03000 

IF(J.EO.N)  GO  TO  87 

03100 

UCF=1 . 0/J.iX>f:)i<2- 1 . 0/DT-FO .  5)?:D  ( J ) 

03200 

ACF  =  j  .0/DX)S:*2+l  .0/DT  +  0.5«D(  J) 

03300 

CTG=T(2?J+1  )  +  T(2?  J~1 )  +  (  ( 1  r  J+1  )  +  T(  1  ?  J--1 ) 

03400 

TN  (  J ) 0 .  GjUCTG/  ( DXS.)f:2«ACF )  -UCF^T  ( 1  ?  J )  /ACF 

03500 

ROT=TN< J)-) (2r J) 

03600 

AfCOT:  ADS  (ROT) 

03700 

IF (AROT. I  F. 0.000001)  GO  TO  77 

03800 

NT  =  1 

03820 

77 

IF( J.NF.Nl )  60  TO  80 

03840 

DNP=0.25*D( J) 

03860 

TlP=0.25#T<lr J) 

03900 

80 

CONTINUE 

03920 

87 

DNP=riNP+0,75«).i(N) 

03940 

TlP=T1P+0.75«r (IrN) 

03960 

TMl =0 . 5« ( T ( 1 f  N1 )  +  T ( 2 »  N1 ) ) 

03980 

TUP=  (  TMl -0 . 5«T  ( 1  r  N )  )  /DX-0 . 5«J.iX)Jc  ( 0 . 25*1  ( 2  ?  N1 )  --T1 P )  /DT 

04000 

1 

-CC*  (  0 . 5*);i  ( N )  *  (  ( 1  r  N )  +  J.iX*DNP*  ( 0 . 25*11-11  +0 . 375*T  ( 1 » N ) )  /B ) 

04020 

Tl.P=0 . 5/DX+O.  3/5*riX/DT  +  CC*  (  0 . 5*1)  ( N )  +0. 375*DX*DNP/B) 

04040 

TN<N)=-TIIP/»LP 

04200 

DO  110  J=1?N 

04300 

T(2?  J)=-TN<  J) 

04400 

1 1 0 

CONTINUE 

04450 

IF(NT.Ff?.0)  60  TO  145 

04500 

NIT=:-NIT  +  1 

04600 

GO  TO  58 

04700 

145 

IK=IK+1 

04750 

TIME' TlMh+DT 

04800 

WRITF<1.»53>  IK»NIT»riMh  »  <  f(2r  J)f  J-1»N) 

r:ar-  rr  .  i  n  mm» 


T53 

T(l, J)"T (2» J) 
CONTI  K'UE 

IF < IK. LI. M)  60  TO 

STOP 

END 


05000 

05100  1 

05200 
05300 
05400 


APPENDIX  H 


HAY 

Tuesday^  November  13)'  1979  21J32J51 


(?TY  SPLI8 

♦  1 

FOR 

00100 

C 

SPII8  — HKAT  TRANSFER  WITH  RADIATION, AN  UNSTEADY, 

00200 

c 

ONE  DIHENSIONA!  FIN  PRODI  FM ,  NOUh HBifR  13,1979 

00300 

c 

ERROR  METHOD 

00400 

DIMENSION  T(2,101)rD(.101  )  ,  TN ( 101 )  ,  AT  ( 101 )  ,  DT  101  ^  ,  DEL  ( 10 

1) 

00500 

OPEN  ( 1JN.TT=:  .1. ,  A(:OE  SS^'  '  BE  01  N  '  ,  Din!  OOUE  ) 

00600 

READ  (1,10)  GC ,  Cr, ,  K(; ,  TR ,  P 

00700 

RE AD (1,20)  M,N 

00800 

10 

FORMA! (8F10.5) 

00900 

20 

FORMAT  (4.1 10) 

01000 

Nl-N-1 

01100 

DX=0C/FL0AT(N1) 

01200 

DT  =  D>;«nX/2.5 

01300 

WR I TE  ( 5 , 30 )  GC ,  CC ,  RC ,  TR ,  B ,  DX ,  J.iT ,  M ,  N 

01400 

30 

FORMAT  ( 1  OX ,  7F 1 0 ,  ,  2 1 1 0/ ) 

01500 

DO  50  1=1,2 

01600 

DO  40  J=1,N 

01700 

T(  I  ,  J)::  1 .0 

01800 

40 

CONTINUE 

01900 

50 

CONTINUE 

02000 

IK  =  1 

02050 

IP=1 

02100 

NIT=1 

02150 

TIME=  0.0 

02200 

WRITE(5,53)  IK, NIT, TIME  ,  (Td,  J)  ,  J- 1  ,N) 

02300 

53 

FORMAT ( 2X , 2T5 , FH . 5 , IX , 5 ( E  9 . 5 , IX ) / , 20 ( 21X , 5 ( F9 . 5 , IX ) / ) / ) 

02400 

56 

NIT  =  1 

02500 

58 

NT=0 

02600 

AT(2)=1 ,0 

02700 

BT(2)= 0.0 

02800 

DO  80  J:=2,N 

02900 

TM=:0 . 5«  (  T ( 1  ,  J)  +  T(2,  J)  ) 

03000 

IK  J )  =  1 . 0+RC>f;  (  rM««3+4  .  ()«TRSTM**2  E6 . 0*TR**2*TMF4 . 0»TR)K)K3 ) 

03100 

CD  J=RC«  ( 1  <  5«TM««2+4 . 0)K  !  R«TH+3 . 0^;TR>i;:J;2 ) 

03200 

IF(J.FR.N)  GO  TO  87 

03300 

UCF  =  1  .0/DX!f:>!<2-l  .0/DT+0.5*D(  J) 

03400 

ACF=1  .o/DX>!;>f;2+i  .o/Dr+o.5^;D(  J) 

03500 

CT&  =  T(2,  J+1  )  +  T(2,  J~1  )+T(  1 ,  JEl  )-!T(  1  y  J--1 ) 

03600 

TN  ( J )  •■=  0 . 5)!;CTG/  ( DX««2)f:Ar;F )  -UCFJf.T  ( 1 ,  J )  /ACF 

03700 

ROT=TN( J)-r(2, J) 

03800 

AROT)=ABS(RnT) 

03900 

CA=2 . 0«DX««2«:ACF  +  0 . 5^:  ( DI  >;:CTG+4 . 0;4:DX**2 )  )KCDJ 

04000 

1 

/(DT^ACF) 

04100 

CB=2.0>l:DX««2)f!ACF 

04200 

AT(  J+1  )::CA-1 .0/A  l  (  J) 

04300 

BT  ( J+ 1  )  :-BT  (  J )  /AT  (  J )  -(;BH<RnT 

04400 

IF (AROT. IF. 0.000001 )  GO  TO  77 

04500 

NT=1 

04600 

77 

IF(J.NE.Nl)  GO  TO  80 

04700 

DNP=0 , 25«).i  (  J ) 

04800 

CDN1=^  CDJ 

04900 

T1NP=  0.25«T(1 , J) 

05000 

80 

CONTINUE 

05100 

87 

if(nt.e;o.o)  go  to  145 

05200 

IiNP=DNP+0.75#D(N) 

05220 

T 1 NP'  T .1  NP+0 . 75«T  ( .1  ,  N  > 

05240 

TMl =0 . 5« ( T ( 1 , N1 ) +T ( 2 , N1  )  ) 

05260 

TUN=(  IM1-0.5*T(  1  ,N)  ) /DX~0 .  EDftDX*  (  0 . 25*!  ( 2 ,  N1 ) -TINP  )/DT 

05280 

mtm 

1 

-CC*(0.rt*r.KN))!t  1  ( 1  ,N)+DX«JiNP*(0.25*TMH0.37G*T(  1  ,N)  )/D) 

/■ 


05290 

TLN=0  ?5/)m0 . 3  753ScrT>:/J.»T+C(;«  ( 0 . 5«ri  ( N )  +0 . 375»DNP#DX/B  > 

05295 

R0T=-TllN/ri..K-T<2»N) 

05300 

SGRU=o«  5/DX-o « .1 2\:mix/ yyr-0 , ( o ,  3753(fT  (i,n)  *cdni 

05400 

/B 

1 

+0 . 5)S:DKK  +  0 . 12r;«  (  f  <  2  r  K'l )  +•(  ( 1  r  K'l  >  )  ^CnK'J  +0 . 375*T  ( 2 » N )  *CDN1 ) 

05500 

SGKEI--  0 . 5/);iX+0 . 375«criX/):n  +nc«  ( 0 . 5«1I  < N  )  +0 . 375«:riX)*;riNt--/B 

05600 

1 

+  ( 0 . 5+0 .  ?8125)i<nX/B )  «r  ( 2 r  K' ) )5;(;riJ+  ( 0 . 5*7  ( j  r  K )+0 ♦  2B125* 

05700 

2 

DX*T  ( .1  r  K* )  /B+0 . 09375«J.IX)};  (  T  ( 2  r  K'3  >  +  (  ( 1  r  K'l )  )  /B )  «Gn J ) 

05800 

SGR=SGRU/SGRIi 

05840 

PRL  =  0 . 5/DX+0 . 375«r.iX/).iT  +  CC*  ( 0 . 5!*!D  ( N )  +0 , 375*DNP*DX/B ) 

05900 

DEL ( K'l  )  - ( -BTC K' >  +ROT!«;PR(.  /SGRB )  / < -SGR+ AT ( N ) ) 

06000 

DEL  ( ?%' )  =  AT  ( K' )  !f:riEI.  ( R1  )  +BT  ( N ) 

06100 

N2=N1-1 

06200 

DO  93  J--=2,N2 

06300 

JR=N~J 

06400 

JR1=N+.1-J 

06500 

DEL(JR)^  (DEI  (JK1)-PT( JRl) )/AT(JRl) 

06600 

93 

CONTINUE 

06700 

DO  110  J=2»N 

06800 

T(2»  J)--T(2r  J)  +  J.iEI..(  J) 

06900 

110 

CONTINUE 

07000 

NIT=NIT+1 

07100 

GO  TO  58 

07200 

145 

IK=IK+1 

07220 

TIME=T1ME+DT 

07240 

IP=IP+1 

07260 

IF(IP.NE.ll)  GO  TO  149 

07280 

IP=1 

07300 

URITE(r:.r53)  IK  >  NI T  r  TIMH  r  (  7  ( 2  r  J )  r  J-^1 » N ) 

074  00 

149 

no  153  J=1»N 

07500 

T(l» J)'  T(2f J) 

07600 

153 

CONTINUE 

07700 

IFdKJ.T.H)  GO  TO  56 

07800 

STOP 

07900 

END 

& 


LOGO 

Killed  Job  38>  User  FMSHEN»  Account  f  TiY  36, 

3t  13-NOV-79  21:3'1:0Aj  Used  0:i2:23  in  A:11;50 


